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1 Introduction 



Localization has been recognized as one of the most powerful techniques in analyzing supersymmetric 
quantum field theories which allows one to calculate some quantities exactly. The localization method of 
[T] was applied to A/" > 3 Chern-Simons-matter theories on in [5]. This result, as well as its A/" = 2 
generalization [3] [4| , can be used to calculate some physical quantities exactly. A recent example of such 
quantity can be found in [5] [6] . 

In many cases, the resulting localization formula for, say, the partition function is still complicated, 
even though it is given as a well-defined finite-dimensional ordinary integral. One way to extract some 
information from the formula is to take the planar limit of the field theory under investigation by which 
the integral can be evaluated exactly in terms of the saddle-point. The analysis of the saddle-point can 
then be performed by the traditional matrix model techniques. The analysis for ABJM theory [7] was 
performed and the planar resolvent was obtained explicitly in [5]. The technique used in [5] was based 
on a relation of the localization formula for the partition function to a topological string theory, but this 
does not seem to be easily generalized to other Chern-Simons-matter theories. A simpler technique was 
introduced in [9] (see also [10] for earlier discussion) which could extract the leading behavior of the free 
energy and the vev of BPS Wilson loops in the large N limit. This technique turned out to be easily 
generalized to various theories including J\f = 2 ones [H] [12] [13] [M] [15] , but a systematic analysis of the 
1/A^ corrections would not seem to be straightforward. 

Later, it was found in [TB] that all the 1/N corrections in the partition function of ABJM theory can 
be explicitly summed up, resulting in a simple function. This surprising result was rederived in |17j by 
a completely different technique, that is, by using a relation of the localized partition function to the 
partition function of a particular free Fermi gas. This Fermi gas picture was also applied to more general 
theories in [T7][T5]P^. 

On the other hand, the traditional analysis of the saddle-point equations in terms of the resolvent 
has been developed in [10] [S] [12 ; but the range of applicability has been rather limited. However, even 
though the progress based the Fermi gas picture is so remarkable, it seems that there are still some 
reasons to investigate the planar resolvent. One possible advantage would be that the planar resolvent 
can provide results which arc exact in the 't Hooft coupling. This means that, if a Chern-Simons-matter 
theory has a string theory dual, the planar resolvent provides the complete information of the worldsheet 
theory with genus zero. In the Fermi gas picture, it is rather difficult to obtain some non-perturbative 
corrections in the sense of the worldsheet theory. Sooner or later, the Fermi gas picture would be improved 
to overcome such difficulties. For a while, the traditional method could be a complementary method to 
analyze Chern-Simons-matter theories. 

In this paper, we investigate J\f = 3 \J{N)k Chcrn-Simons theory coupled to an arbitrary number of 
adjoint matters. We call this family of theories as A/" = 3 Chern-Simons-adjoint theories. An J\f ~ 2 
version of such a theory was recently discussed in [18] and some qualitative properties of the theories were 
found. Our theories have larger supersymmetry, and the analysis of the theories is easier. We determine 
the planar resolvent of the theories. The resolvent is given as an integral of a combination of elliptic theta 
functions. This integral expression of the resolvent allows us to analyze the analytic continuation of the 
't Hooft coupling and the vev of a BPS Wilson loop [13]. The eigenvalue distribution obtained from the 
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resolvent has a similar property observed in |18] . 

This paper is organized as follows. In section ^ we explain an ansatz for solving the saddle-point 
equations. This ansatz is applied to A/" = 3 pure Chern-Simons theory in section [3] to check the validity 
of the ansatz. Then, Af = 3 Chern-Simons-adjoint theories are analyzed in section |4l We give integral 
formulas for the planar resolvent, the 't Hooft coupling and the vev of the Wilson loop as functions of the 
position of the eigenvalue distribution. Section [5] is devoted to discussion. There are several appendices. 
Appendix El reviews a solution of a homogeneous equation studied in [21] [25] . Appendix |B] contains a 
proof of an identity which is used in our analysis. Appendix [C] discusses the conditions by which the 
solution of the saddle-point equations is determined uniquely. Appendix [D] contains some details of the 
calculations which determine the planar resolvent. 



2 Matrix models 



In this paper, we consider J\f = 3 Chern-Simons-adjoint theory, that is, Af ~ 3 Chern-Simons theory 
with the gauge group V{N) and the Chern-Simons level k coupled to 2n matter fields in the adjoint 
representation of V{N). The number of the adjoint matters is chosen to be even so as to preserve Af = 3 
supcrsymmetry. Some basic facts on A/" = 3 Chern-Simons-matter theories can be found in [23] . 
The localization formula [2] for the partition function oi Af = 3 Chern-Simons-adjoint theory is 



exp 



i=l 



n.<, sinh^ 



(2.1) 



The normalization of Z here is different from, for example, the one in |17| . This is not an important 
modification since the overall constant is irrelevant in the following analysis. 

The vev of a BPS Wilson loop [23] can be also given as a finite-dimensional integral as [2] 



(W) 



1 f ^ 

z y n '^"^ 



i=l 



n,<,sinh2 



cosh^ 



1 ^ 
N ^ 



(2.2) 



2.1 Saddle-point equations 

In the large limit, the integral in the localization formula ()2.1|) can be evaluated exactly by the saddle- 
point configuration {wi} which satisfies 

—u, = ^coth^^-n)Jtanh^^, (2.3) 



3 



where i and j run from 1 to . and —ik is replaced with k, understanding that k can take a complex value. 
These equations are invariant under the sign flip Ui — )■ — w^. Therefore, the saddle-point configuration is 
expected to be symmetric under the sign flip. That is, an equality between two sets 



{ui,--- 



UN } 



{ -Ul, • • • , -UN } 



(2.4) 



is expected to hold. 

When k is real and positive, the equations (|2.3p can be regarded as equations of the balance among 
forces. The left-hand side represents an external linear force acting on i-th eigenvalue Ui which confines 
it around the origin. The right-hand side represents the sum of forces between Ui and the others. The 
first term in the right-hand side provides a repulsive force, while the second term provides an attractive 
force. It is natural to expect that the values of Ui are distributed around u = 0, at least when k is large 
enough so that the external Gaussian potential is dominant. 

The nature of the interaction among the eigenvalues depends of n. The case n — corresponds 
to pure Chern-Simons theory whose large N solution is known [26|[27]. There are long-range repulsive 
interactions, and due to them, the eigenvalue distribution can be extended largely when the 't Hooft 
coupling A :~ N/k is large. For the case n — 1, there is no long-range interactions, like in the case of 
ABJM theory. The large N solution for this case was obtained in [22]. In this case, the width of the 
eigenvalue distribution can be large for large A, but its scaling with respect to A is different from neither 
the one in pure Chern-Simons theory nor in ABJM theory. For the remaining cases n > 2, there are 
long-range attractive interactions. It it expected that the eigenvalues could not be largely extended in 
these cases. 

There is a quick way to guess whether the support of the eigenvalue distribution can be large if A: > 0. 
If Ui is largely separated from most of the other eigenvalues, the hyperbolic functions in (|2.3p may be 
approximated by the sign function |10j . The resulting equations 



would make sense if n < 1, while they look inconsistent if n > 1. This suggests that the eigenvalues could 
be broadly distributed in the former case, but the eigenvalues would be always confined in a finite range 
in the latter cases. 

The level k in (j2.3p must be replaced with —ik with an integer k so as to obtain some results relevant 
to A/" = 3 Chern-Simons-adjoint theories. This will be achieved via an analytic continuation of k after 
the solution of (|2.3[) for a positive real k will be obtained. 

It is convenient to introduce new variables Zi := e"'. In terms of them, the equations (|2.3p can be 
written as 




(2.5) 




(2.6) 



where t := 27rA. The symmetry (|2.4p implies 



{zi,--- ,zn} 



-1 
1 ' 



(2.7) 



The eigenvalues {zi\ are distributed around z = 1. 
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A traditional strategy to solve (|2.6p is to define a suitable complex function usually called the resolvent. 
A convenient choice of the resolvent in this case is 

i 

For large z, v{z) is expanded as 

v{z) = t + 2t{W)z-^ + 0{z-^), (2.9) 
where (W) is the vev of the Wilson loop (|2.2p . In the large N limit, it is given as 



1 ^ 

{W) = ^E^- (2.10) 



The symmetry property (j2.7p implies 



= -^^(2), (2.11) 
where 7^ is assumed. This identity then implies that 

w(0) = -t (2.12) 

which is compatible with the definition (|2.8p . 

As usual in the matrix model literature, we assume that the large N limit would be described by a 
continuous function p(x) for the eigenvalue distribution. The meaning of p{x) is 

1E^(^--^^) ^(^) (2.13) 

i 

in a suitable sense. The integral of p{x) is normalized to be 1. According to the discussion above, the 
support of p{x) consists of a single segment [a, b] with 

ab = 1, (2.14) 

which comes from the symmetry (|2.7p . The resolvent (|2.8p is written as 

v(z) = t / dxp{x)^-^. (2.15) 
J z - a; 

In terms of v{z), the equations (j2.6p can be written as 

2 logy = v(y + iO) + u(?/ - iO) - 2nij(-?/). (2.16) 
Here and in the following, the variable y always takes a value in [a, 6]. 

The expression p.lSp implies that v{z) has the following analytic properties. 

• v{z) is holomorphic on C\[a, 5]. 

• v{z) is finite at z = a,b since p{x) vanishes at a: = a, b. 

• v{z) is finite at infinity. 

It is shown in Appendix ICl that the solution of (|2.16|) with these analytic properties is unique. 
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2.2 Ansatz for the resolvent 



To solve the equation (|2.16p . we make the following ansatz for v{z). We assume that v{z) can be written 
as ^ 

viz) = I d^viz,0, (2.17) 

where v{z^£^) is assumed to be a solution of 

= v{y + iO,0+v{y-iO,0-'2nv{~y,0- (2-18) 
We further assume that v{z,£_) has the following analytic properties. 

• w(2,^) is holomorphic on C\[a,b]. 

• v{z, ^) is finite at z = a, b. 

• is finite at infinity. 

Note that a,b are common for w(z,^) with any — oo < ^ < 0. The argument in Appendix [Cl shows that 
such av{z,£,) is also unique. 

It is easy to show that v{z) given as ()2.17p satisfies ()2.16p . Moreover, as long as v{oo) turns out to 
be finite, p.l7p has all the required analytic properties for v{z). Therefore, due to the uniqueness of the 
solution of (|2.16p . it turns out that one can obtain the solution of (|2.16p by solving (|2.18p . The equation 
()2.18|) is a bit simpler than p.l6p since the left-hand side of ()2.16|) has a branch cut while that of (|2.18p 
has only a pole. 



3 J\f = 3 Pure Chern-Simons theory 



The simplest case n — corresponds to M ~ 3 pure Chern-Simons theory. The planar solution for this 
theory is known |26] [27j . In this section, we show that our ansatz ()2.17|) can reproduce the known results 
for this theory. 

3.1 Planar resolvent 



For the case n = 0, the equation (|2.18|) is simply 

2 y-l 



Define a function w(z,^) such that v{z,^) is written as 



v{y + iO,0 + viy-iO,0- (3.1) 



v{z,0 = -j^^+u;iz,0- (3.2) 
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Then, p.ip implies that uj{z,£,) satisfies 

= w(y + ^0,O+^(y-^0,e). (3.3) 
The general solution of this equation is of the form 



ioiz,0 = A{z,^)^{z~a){z-b), (3.4) 

where A{z, ^) is a meroniorphic function on C. Note that the finitcness condition at z = a,b was used. 

Since v{z,^) is assumed to have no poles, cj(z,^) must have a pole at z = ^ with the residue 1. This 
implies that A{z,^) must have a pole at z = ^. Therefore, A{z,£^) is of the form 

AizA) = (3.5) 

where Ai{z,^) is an entire function. Requiring the finiteness of v{z,^) at infinty, Ai{z,£,) must be finite 
at infinity, implying that Ai{z,^) is independent of z. The ^-dependence of Ai(z,^) is then determined 
by the required residue of the pole of a;(z, ^) at z = ^. It turns out to be 

M^,0 - \., (3.6) 

It is concluded that u(z,^) can be written as 



t 1 z-l 1 ^{z-a){z-b) 

"'^•f ^ -——(^— vi(-^}((-»y 

The resolvent v(z) is obtained as the integral (|2.17p of v(z,£^). 

The identity (|2.1ip expected from the general argument can be verified explicitly. Using ah = 1, one 
finds 



1 1-z 1 ^(1 - az)(l - bz) 



= -rMz,r'), (3.8) 

implying v{z) = —v(z~^). 

Let us check whether the solution p.7p is the correct one. There is the following integral expression 
of the resolvent 1221. 



^ 2 / ^ log(^) Viz~a)i z_± 
Ja z-x ^\{x-a){x-b)\ 

Note that the definition of the resolvent in [H] is slightly different from (|2.15p . Introducing a new variable 

(3.10) 



{x - l)r; + 1 ' 



one can rewrite the integral in p.9p as 



dx log x \/ iz — a)(z — b) ^ I '^^ x — I 1 \/ {z — a)(z — b) 

2^z-x ^x~a){x-b)\ ~ Jo Jc^^H^- + lz-x^(x- a){x ^ 

^ dT] (f ^^-i ^'v \/ {zrj ~ a,i)[zn - br,) ^^^^^ 



7 



where C and Cn are contours encircling the segment [a, h] and [a,,, 6^], respectively. Note that the following 
identity was used 



1 z-l 



logz = / d,J^^^,- (3.12) 







The ic^-integral provides the resolvent of the Gaussian matrix model with the branch points and 
Using the explicit form of the resolvent of the Gaussian matrix model, one finds 



2 ■ dx \ogx ^{z - a){z - h) ^ /"^ ^ 



(3.13) 



By the change of integration variable = (1 — ^ , one can show that the integrand coincides with p.7p . 
This confirms the validity of our ansatz ()2.17p for the resolvent. 

Interestingly, the integral form p.l3p suggests that there is another way to obtain the solution of 
(|2.16p for the case n = 0. Using the identity p.l2p . one finds that the solution of (|2.16l) can be obtained 
from the solution of 

2y„ = w(^/ + ^0,^?) + ^;(^/-^0,?7). (3.14) 

Define v{z,r]) such that 

v{z^,T]) = v{z,r]). (3.15) 

Then, v{z,ri) satisfies the saddle-point equation for the Gaussian matrix model with the branch cut on 
[a,,, 6,,]. Due to the required analytic properties of v{z,r]), v{z,rf) must be holomorphic on C\[a^,6^]. In 
addition, v{z,t]) must be finite at infinite since 

1 — 7? 

z„ -> 00 ^ z ^ -. (3.16) 

1] 

These conditions imply that the resolvent is given as in (|3.13|) . 

Unfortunately, this method cannot be applied to the general case n > 0. If one writes down the 
saddle-point equation following the above strategy, one obtains 

2y>7 = v{y,, + i0,i]) + v{y^-i0,r])-2nv{{~y)r,,r]). (3.17) 

Since {—y),j 7^ ~2/r;, one cannot simply borrow the results in, for example, [^[25]. 



3.2 Analytic continuation 



The explicit formula p.7p for f(z,^) and the relation t = — w(0) provides an integral formula for the 
't Hooft coupling 



tia) 



-00 C 



1 



1 



(3.18) 



as a function of a. This integral is well-defined for < a < 1. Note that \/{^ ~ a,){£, — b) is defined such 
that 

Vi^-aM-b) ^ ^, (e ^ -foo) (3.19) 
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Figure 1: The deformed contour for the 2-k rotation of the phase of a (accompanied by — 27r rotation 
of the phase of b) (left). The contour on the right is an equivalent contour. The broken lines lie on the 
second sheet on which the integrand has a pole at the origin. The branch cut is shown by blue lines and 
the simple pole at ^ = 1 is shown by the crosses. 

so that the integrand is finite at ^ = 0. The finiteness p.l8|) confirms the finiteness of v{z) at infinity 
which is necessary for the uniqueness of the solution of (|2.16p . 

It is apparent from (|3.18p that t{a) is always real for < a < 1 and t{l) = 0. An arbitrarily large 
t(a) may be realized if there is a limit of a such that the integral diverges. Indeed, a — > is such a limit, 
and the divergent contribution comes from the integration near ^ = 0. One finds, for small a, 

t{a) = -loga + 0(l). (3.20) 

To obtain some results relevant to pure Chern-Simons theory, an analytic continuation must be 
performed so that t(a) becomes purely imaginary. This can be done by regarding a as a complex variable, 
and then extend the domain of t(a). Since the condition a6 = 1 is necessary for the integral p. 181) to be 
convergent, the analytic continuation must be done while keeping this relation. Without loss of generality, 
one may assume \a\ < \b\. Therefore, t{a) can be regarded as a function defined on a region 

:= {a e C I < |a| < l,a 7^ -1}. (3.21) 

A point a = — 1 was eliminated since for this value the two branch points = a,b are merged together 
on the integration contour, resulting in the divergence of t{a). Since the integration contour is pinched 
by the branch points, there is no way to deform the contour to avoid the divergence. 

The function t{a) is in fact a multi-valued function on D. Indeed, this can be anticipated from the 
log-behavior ()3.20|) . It turns out that the integral expression (|3.18|) is very useful to analyze the analytic 
structure of t{a). Consider the phase rotation of a. The branch points of the integrand move around 
in the complex ^-plane, and due to this, the integration contour must be deformed to avoid the branch 
points. After the 27r rotation, the integration contour becomes the one depicted in Figure [T] Figure [1] 
also shows an equivalent contour. Note that the integrand is smooth at ^ = on the first sheet, but has 
a simple pole with the residue 2 on the second sheet. Therefore, one obtains 

t{e'^''^a) = t{a) - 27ri. (3.22) 

This is consistent with the log-behavior p.20p . 
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Figure 2: The unit disc in the a-plane (left) is mapped into an infinite strip with the width 2-k in the 
t-plane (right). The red curves corresponds to \a\ = i. 

As a map, t{a) maps D into an infinite strip in the complex i-plane, as depicted in Figure [21 The 
only way to obtain a large imaginary part of t{a) is to rotate the phase of a many times. The modulus 
\a\ should be chosen to be of order 0(1) so as to make the real part of t(a) vanish. As a result, the 
eigenvalue distribution, which is given by the branch cut of v{z), has a length of order 0(1) even when 
Im(t(a)) is chosen to be large. 

All the results obtained above can be easily derived from the explicit relation [22] 

v^+^ = 2e^*. (3.23) 
It will turn out in subsection 14.31 that a similar analysis can be performed for the general cases n > 2. 



3.3 Wilson loop 



As shown in subsection 12. 1[ the vev (W) of the Wilson loop can be derived easily from the resolvent. 
One obtains an integral representation 



2t{W) 



a+b 
2 



1 



(3.24) 



This integral is well-defined. Since the integrand has no pole on the whole Riemann sheet, the integral 
does not change by the 2ti rotation of the phase of a. In other words, the integral defines a single- valued 
function on D. The integral diverges in the limit a — > as 



2t{W) = 2a-^ +0{l). 



(3.25) 
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On the other hand, the integral is finite at a = —1. 

The vev (W) itself is a multi-valued funetion on D sinee t is mufti- valued. For a fixed \a\, \{W)\ 
decreases as the phase of a decreases, or Im(i) increases, since the integral p.24|) is bounded on \a\ = 
const. 

One can check the validity of the above results by examining the explicit formula 

2t{W) = J^r'^^T,, - (3.26) 



Note that, using p.23p . one finds 
which is the large N limit of the expression obtained in [2] . 



{W) = e^*^ (3.27) 



4 J\f = 3 Chern- Simons-adjoint theories 

Let us return to the general case. The equation to be solved is 

2 y-1 



v{y + iO, + v{y ~ iO, i) ~ In vi-y, 0- (4-1) 



In the following n 7^ 1 is assumed. The case n = 1 was studied in [35] based on a different technique. 
Let be defined as 

As a function of z, f{z,£,) has a simple pole at 2 = ^ with the residue 1. Then (|4.ip becomes 

-2/(y,0 = viy + iO,0 + viy-iO,0-^nvi-y,0- (4-3) 
Define cj(z,^) such that v{z,^) is written as 

= (4.4) 
— 1 

The first term in the right-hand side is chosen such that 

^fiy,0+nfi~y,0 _^Ji-y.O + nfiy,0 ^ ^/(j;,^. (4.5) 
n — I — 1 

Therefore, uj{z,^) satisfies the following homogeneous equation 

= u;{y + iO,£_) + uj{y - iO,^) ~ 2nuj{~y,^). (4.6) 

Fortunately, this equation was studied intensively in [HjlH] in the context of O(A^) matrix model 
|28j |29] ■ In particular, an explicit solution in terms of the elliptic theta functions was found in [25]. 
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4.1 General solution of the homogeneous equation 



We review the method for solving 

= uj{y + i0)+uj{y-i0)-2nuj{-y) (4.7) 

which was shown in [25) . It is convenient to parametrize n by as 

n = — cos(7ri'). (4.8) 

Since our main interest is in the cases n > 2, the parameter v will be analytically continued to a complex 
value. Specifically, v will be taken to be 

V = l + ii>, i> > 0. (4.9) 

The equation (|4.7p can be written as 

uj±{y + iO) = -e=^"''w^(y- iO), (4.10) 

where uj±{z) are defined as 

u;+iz) := ^' c._(z) := uj+{-z). (4.11) 

2 sm(7rz^) 

There exists a pair of functions G±{z) which satisfy ()4.10p [21] [25]. The properties of G±{z) are 

summarized in Appendix |X| The general solution of (|4.10p is obtained in terms of G±{z). Let S±{z) be 
defined as 

^^^'^ G^y ^'-''^ 

Then, S±{z) satisfy 

S±iy + iO) - S^{y~iO). (4.13) 

The implications of these equations become transparent when they are written in terms of the following 
functions, 

S{z) := 5+(z) + 5_(z) = S+{z) + S+{-z), (4.14) 

S{z) ;= S+{z)-S-{z) = S+{z)- S+{-z). (4.15) 

Note that S{z) is an even function while S{z) is odd. They satisfy 

S{y + iO) = S{y~iO), S{y + iO) = -S{y-iO). (4.16) 

These equations imply 



Siz) = 2A{z^), Siz) = 2zB{z^)^{z^ -a^){z^ -b'^), (4.17) 

where A{z'^) and B(z'^) are meromorphic functions on C. In terms of these functions, u}±{z) can be 
written as 



w±(z) = ^(z^) ± zS(z2)^(z2 _ a2)(^2 _ 52) G±(z). (4.18) 
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Since uj-{z) is obtained as lo+{—z), it is enough to determine uj+{z). 



As reviewed in Appendix lAl G-^{z) has singularities at z = a, 6 where G-j-{z) behaves as (z — a) 2 
and {z — b)~^ , respectively. Assuming that w+(z) is finite at these points, it must be of the form 



uj+{z) = Ai(z2)^(z2 - a^){z^ - 62) + zBiz^) ^{z^ - a^){z^ - b-^)G+iz) 



(4.19) 



G+{z) has a simple zero at z = e (e is determined by a and b. See Appendix This allows both 
Ai{z'^) and B{z^) to have simple poles at z = e even when w+(z) is assumed to be regular at z = e. 
Then, one can write w+(z) as 



cj+(z) = 



,, V(z2-a2)(z2-62) 

^2(2 ) :5 — :5 + — -5 



\/(z2-a2)(22_52)c<^(^)^ (4 20) 



Let us assume further that a;+(z) is regular at z = — e. Then, the pole at z = — e in the above expression 
must be canceled. Rewriting Bi{z^) as 



i?i(z2) = lA2{z^) + B2{z^), e := ^ [e^ - a^){e^ - b^) , 
e 



one obtains 



where 



uj+{z) 



A2{z^)g{z) + z 



B^iz^V 



^{z^-a^){z^-b^)G+{z), 



v/(z2-a2)(z2-62) + fz 



(4.21) 



(4.22) 



(4.23) 



The regularity at z = — e then implies B2{z'^) = (z2 — e^)B^{z^). 



In summary, it has been shown that the solution of the equation (j4.10p which is finite at z = a, 6 and 
regular at z ~ ±e has the general form 



Lj+{z) = [A{z^)9{z) + zB(z^)] V(^2 _ a2)(z2 - 62)G+(z), 



(4.24) 



where A{z'^) and B{z'^) are meromorphic functions regular aX z — a,b,zte. 



4.2 Planar resolvent 



The function uj{z,^) defined in (|4.4p depends on the parameter ^. Recall that the parameters a,b are 
assumed to be independent of ^. Therefore, the solution in the previous section can be used for uj-f.{z,^) 
simply by allowing the functions A{z'^) and B(z'^) to depend on ^. That is, aj+(z,^) can be written as 

w+(z,0 = [A{z\ Oaiz) + zBiz\ 0] Vi^^ - - &^)G+(z). (4.25) 

It will be shown below that A(z2,^) and i3(z2,^) are determined by the analytic properties of aj+(z,^) 
expected from those of ?;(z,^) assumed in subsection 12.21 
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The definition (|4.4p oiuj{z,^) implies that is related to v{z,^) as 

= v+{z,0 - ^^J_ ^) (e-^"V(^,C) - /(-z,g)) , (4.26) 

where is defined from t;(z,^) as in (|4.1ip . This relation then implies that cj+(z,^) is finite at 

infinity. This also implies that must have simple poles at z = ±^ with the residues 

Res^ c.H.(z,0 = - 2(^2'„i) «~^"'^ (4-27) 
Rcs„5C.+ (z,0 = - 2(„2'„i) ^""^- (4-28) 



These poles must come from A{z'^,^) and Therefore, they can be written as 

(4.29) 

where v4i(z^,^) and -Bi(z^,^) are entire functions of z. The behavior of near infinity is then 

w+(z,e) ^ »z-Mi(z2,a+zBi(z2,0, (4.30) 

where 

G+(z) 1 + 0(^-2) (4 3^) 

z 

was used. See Appendix El The finiteness of uj^{z,£^) at infinity implies 

M^^O = c(e), = diO- (4.32) 

The constants c(^) and d{^) are then determined by requiring that aj_|_(z,^) has the correct residues 
(|4:27))(|4:28)) at z = One finds 

ciO = rf(0 = ^(;^r^G(r^), (4.33) 

where the following identity was used, 

eg{z)G+{z) = z'^G+{z-^), (4.34) 
which is valid when a& = 1 is satisfied. See Appendix [B] and [Dl for the details. 



Now, aj+(z,^) has been determined completely. It can be written in the following form 

1 

2(n2 - l)e 



^+{^,0 = / Uz,o-r'^iz-\r') , (4.35) 



where 



-(z,0 V(^^ ^^"■^)(^^^ ^\ GiOz-'G^{z~'). (4.36) 



z2-e 

It is straightforward to verify the identity ()2.1ip by noticing 



fiz-\o = -rV(^,r'), (4.37) 
io+{z-\o = -rV(2,r'). (4.38) 

14 




Figure 3: The contour in Figure [T] is drawn in the it-plane (left). The upper and the lower edges of the 
reetangles are identified. The crosses correspond to ^ = 1. An equivalent contour is also drawn (right). 



One might be worried that the results obtained above would depend on a particular choice of v. 
Indeed, since n = 1 is excluded in our calculation, there are at least two ways to continue n from 0, which 
was studied in section |3l to the cases n>2. Understanding (|4.8|) as a map from the complex n-plane to 
the complex i/-plane, the upper-half n-plane is mapped into 

R+ := {v e<C\Q< Rc{u) < 1, Im(i/) > }, (4.39) 

and the lower-half n-plane is mapped into 

R- := { z/ e C I < Re(i/) < 1, Im(i/) < }. (4.40) 

These two choices correspond to taking i> to be positive, as in (|4.9p . or to be negative. It can be checked 
that both prescriptions for the analytic continuation provide the same a;(z) by noticing that the only 
change due to the change of the prescription is the sign flip of G{z) . The independence of the prescription 
is also a consequence of the uniqueness of the solution of ()4.7|) based on an argument similar to Appendix 



4.3 Analytic continuation 



Recall that the 't Hooft coupling t is obtained as a function of a through the resolvent as 

tia) = 



MO) 
rO 



c^(o,e) 



1 1 



The explicit expression of (|4.35|) gives 



t{a) = 



1 



7^-1 J-oo £. 



2e sin ^ 



n-ie(e-l) 



-GiO 



(4.41) 



(4.42) 



Note that this integral expression is quite similar to (j3.18p for pure Chern-Simons theory. 
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Figure 4: The red curve is the plot of et{a) against a, and the blue curve is the plot of e against a. 
n ~ 5 is chosen. 



This integral is well-defined for < a < 1- Indeed, for large ^, the integrand behaves as ^ ^, while at 
^ = the function in the square brackets vanishes since 

G(0) = 2sin— G+(0) = 2iesin— . (4.43) 

The finiteness also implies that v{oo) is finite, as expected. 

The integral ()4.42|) defines a function on D. As in the case of pure Chern-Simons theory, t{a) is a 
multi- valued function. One may expect this property from the fact that G{z) has the branch cut on [a, 6], 
the same position as ((z — a){z — h))~^ . Due to this analytic structure, the effect of the phase rotation 
of a by 27r is determined by the deformation of the integration contour which is exactly the same one 
depicted in Figure [1] It is more convenient to consider the contour in the u-plane. The deformed contour 
in the u-plane is shown in Figure |3l Note that the vertical contour in Figure |3] is equivalent to a contour 
encircling the segment [—6, —a] in the ^-plane. Since the integrand has no branch cut on [—6, —a], the 
contribution from this contour vanishes. Now, the contour deformation amounts to adding residues of 
two poles, C = 1 E^nd ^ = on the second sheet [u ~ 2K). One finds 



n — 1 



"shT^ 



1 - 1 



In — 1 

= t{a)-2T:i . (4.44) 

71—1 

For the case n = 0, this formula coincides with p. 221) for pure Chern-Simons theory. The singularity at 
n = 1 seems to be compatible with the results in [2^ . 

A plot of et{a) for real a with a fixed n is shown in Figure ID The choice et{a) is due to the fact that 
a factor exists in the integrand which may diverge for some a. Figure |4] shows that t{a) increases 
from as a decreases from 1. The positivity of t{a) for a real a can be regarded as a consistency check 
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Figure 5: t{a) is plotted against n with a chosen such that = 2. t{a) diverges at n = coshyr ~ 11.6 
according to (|4.45p . 



of our resolvent since it is required from the definition. Note that e is also a function of a which may 
vanish even when et{a) is non-zero. This indicates that t(a) diverges to +oo at a = a* where a» satisfies 



See Appendix [XI for the dependence of K,K' on a. In other words, t{a) has a simple pole at a = a*. 
Recalling that e is proportional to the elliptic function STv(yK' ,k) and that K' increases as a decreases, 
one finds that t{a) has an infinite number of simple poles on the real axis in D. 

If a becomes smaller than a,, then t{a) becomes negative. This is certainly not a signal of the 
inconsistency of the resolvent, but the indication that for any real and positive 't Hooft coupling, a is 
always within the range (a*,l]. In other words, the width of the eigenvalue distribution is always of 
order 0(1). This behavior of the eigenvalue distribution fits nicely with the intuition obtained from the 
original saddle-point equations (|2.3p . For the case n > 1, there are long-range attractive interactions 
among eigenvalues, so they would be bound within a finite region. On the other hand, one can check that 
e never vanishes for the case n < \. Due to this fact, a can be arbitrarily close to before t{a) diverges, 
resulting in a largely extended eigenvalue distribution. Correspondingly, there are long-range repulsive 
interactions among eigenvalues for the case n <\. 

It would be interesting to see the dependence of t(a) on n. Figure shows a plot of t{a) as a function 
of n with a fixed a. The value of t{a) for n = reproduces the known exact value 



It turns out that t{a) is positive as n increases from up to some value. Interestingly, t{a) is smooth 



vK' 



2K. 



(4.45) 




(4.46) 
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Figure 6: A small disc centered at a = a* in the a-plane is mapped to a neighborhood of infinity in the 
t-plane. 

at n = 1 although the expression (|4.42l) looks singular. As n increases further, Figure [5] shows that t{a) 
diverges to +oo. This is again due to the presence of the factor e"""^ in (j4.42p . 

The issue of the analytic continuation from a real t{a) to t{a) = 2'Ki\ with a large positive A is rather 
subtle. For pure Chcrn-Simons theory, there is a unique way to obtain such a value of t{a). This is 
because, although t{a) is multi-valued, its inverse function a{t) is single- valued. However, for the case 
n > 1, there are infinitely many values of a for which Im(i(a)) becomes large. One way to obtain a large 
Im(t(a)) is, as for pure Chern-Simons theory, to rotate the phase of a many times to increase Im(t(a)). 
Another way is to choose one pole of i(a), say a ~ a^, and choose one point in the vicinity of a*. Because 
of the pole at a = a,, t{a) maps a small disc centered at a» in the a-plane to a neighborhood of infinity 
in the t-planc. This is depicted in Figure |6l Therefore, by giving a tiny imaginary part to a near a,, 
one obtains a large and positive A. If the latter continuation would be chosen, then observables would be 
insensitive to the analytic continuation as long as the observables are continuous at a = a,. 

The ambiguity for the prescription of the analytic continuation can be also described as follows. The 
function t{a) has many poles on its Riemann sheet, which consists of an infinite number of D glued 
together, and also diverges along the direction corresponding to the infinite phase rotation. Therefore, a 
given large value of Im(i(a)) does not uniquely specify the value of a, and therefore, the planar resolvent. 
One of the most natural choice would be to choose the pole a = a, since it is directly connected to a 1 
or t{a) for which the matrix model is apparently well-defined. 

Note that the width of the eigenvalue distribution is of order 0(1) independent of the prescription of 
the analytic continuation. 
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4.4 Wilson loop 



An integral formula for the vev {W) of the Wilson loop is 



2t{W) 



1 



71-1 



c 



where C is a constant appearing in the large z expansion of G+{z), 

iC 



+ ^+0{z-^). 



I 

z z^ 

The integral is well-defined. Note that the integral decreases fast enough for large ^ since 



2 cos ^+o(r^). 



(4.47) 



(4.48) 



(4.49) 



Since there is no pole in the integrand, the integral defines a single-valued function on D. Note that 
2t{W) may diverge at some a due to the factor e~^, but this divergence is due to the divergence of t, and 
therefore {W) is finite for a non-zero a. 

The large A behavior may depend on the choice of the analytic continuation. If one obtains a large 
imaginary part of t by rotating the phase of a many times, then {W) decreases as 0{t~^) since the integral 
()4.47|) is bounded for a fixed \a\. One the other hand, if one choose an analytic continuation in which a 
is chosen to be a^, plus a suitable infinitesimal imaginary part, then one finds {W) ~ 0{1). 



5 Discussion 



In this paper, wc determined the planar resolvent oiAf — 3 Chern-Simons-adjoint theories. The resolvent 
was given as an integral of a combination of elliptic theta functions. Although the expression looks 
complicated, it is explicit enough to investigate its analytic structure. It was clarified that there are 
various ways to continue the 't Hooft coupling t to a purely imaginary value by which one can hope to 
extract some information on Chern-Simons-adjoint theories. For those analytic continuations, the width 
of the eigenvalue distribution does not grow indefinitely in the large lm{t) limit which can be deduced 
from the saddle-point equations p.3|) . Note that a similar behavior was observed in [18] for an = 2 
version of Chern-Simons-adjoint theories. Accordingly, the vev (W) of a BPS Wilson loop is finite. It 
was observed that different behaviors of (W) may result from different analytic continuations. 

It is interesting to generalize the analysis performed in this paper to more general Chern-Simons- 
matter theories. It is expected that the addition of fundamental matters would be rather straightforward 
since it turned out to be the case for pure Chern-Simons theory [20] and for AB JM theory [21] . Inclusion 
of matters in the symmetric and the anti-symmetric representations would be possible in the following 
way. A pair of the symmetric matters would contribute to the saddle-point equations (|2.3p as the following 
additional terms 

1 -1 + -1 -, \ 

~ 2 / ^ tanh — ^ — tanhui. (5.1) 
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Note that these terms does not spoil the symmetry (|2.4I) . In terms of Zi, these terms can be written as 

1 ^ Z^- Z^ I Z,- Z^^ = 1 ~ ^3 1 ^» ~ 2) 

2^z, + z-' 2z, + z-i 2^z, + z, 2z, + z-i' 

where the symmetry (|2.7|) was used. Therefore, the effects of adding a pair of the symmetric matters 
appear as a shift of n by ^ in addition to a modification of the external force acting on Ui or Zi. The 
effects of anti-symmetric matters is almost the same, but the sign of the extra external force is opposite. 

It would be more interesting to investigate Chern-Simons-matter theories with multiple gauge group 
factors. The analysis in this paper would be able to be generalized to quiver-type theories with two nodes. 
In these theories, one can also consider bi-fundamcntal matters. A systematic analysis of such a family 
of Chern-Simons-matter theories would provide an interesting pattern in the eigenvalue distribution and 
the behavior of the Wilson loop. 

Of course, it would be interesting to ask whether there exists a dual description for those Chern- 
Simons-matter theories. Since the Wilson loop in most of the theories mentioned above would not 
exhibit the exponential behavior, which would be expected if there exists a dual classical worldsheet in 
AdS space, the gravity dual would not be simple even if it exists. The possible dual theory might be 
a higher spin theory, as in [30] [34] [35| , or it might be something else. The detailed study of the planar 
solution may provide some information necessary for this line of research. 

Very recently, it was shown in [31] that the free energy of J\f = 2 super Yang-Mills theory in five 
dimensions scales as , as expected from the conjectured relation to (2,0)-theory in six dimensions 
[32j |33| . Curiously, the same calculation can be performed for A/" = 3 pure Chcrn-Simons theory with 
purely imaginary level, and the free energy scales as . Note that this scaling comes as . In addition, 
it was shown in |36| that a modified partition function has a localization formula whose perturbative part 
coincides with the one for pure Chern-Simons theory. In pure Chern-Simons theory, since the level must 
be an integer, the above leading term becomes the imaginary part of the free energy. This observation 
might shed some light on the physical meaning of the scaling behavior. 
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A G{z) 



This appendix shows some properties of a function G{z) by reviewing the derivation given in |25j . 
The function G{z) is a solution of the homogeneous equation 

G{y + iO) + G{y-iO)-2nG{-y) ^ 0. (A.l) 

It is assumed that G{z) is a holomorphic function on C\[a, 6] and that G{z) has a branch cut on [a,b]. 
Due to this branch cut, the Riemann sheet of G{z) would be non-trivial. To make the analysis of the 
Riemann sheet tractable, it is convenient to introduce a new variable u via the Jacobi's elliptic function, 

z — asn(u, fc), (A. 2) 

where k := —. Note that k in this appendix has nothing to do with the Chern-Simons level. Some points 
b 

on the 2-plane are mapped as follows, 

z = 0, a, b, oo ^ u = 0, K, K + iK', iK' (A.3) 
where K and K' are the complete elliptic integrals 

K = K{k) = t , '^^ K' = k(VT^). 



^(l-x2)(l-/fc2x2)- 

The advantage of introducing u is as follows. The inverse function of (|A.2p is 

dx 



(A.4) 



(A.5) 



/o ^{l-x^){l-k^x^) 

This defines a particular kind of the Schwarz-Christoffel mapping which maps the whole z-plane into a 
rectangular region Rq where Ri for an integer / is defined as 

Ri := {ue€\{2l- l)K < Re(u) < {21 + l)K, -K' < Im(u) < K' }. (A.6) 

Note that the upper and the lower sides of [a, b] and [—6, —a] are mapped on different regions. The 
mapping is depicted in Figure [71 An extension of G{u) :~ G{z(u)) from Rq to the whole u-plane in an 
appropriate way determines G(u) on the whole Riemann sheet. 

The function G{u) satisfies some equations. The first one is the homogeneous equation (jA.l[) . In 
terms of u, this is written as 

G{u + K) + G{-u + K) -2nG{u- K) = 0, (A.7) 

where u € i[0,K']. The second one is the equation 

G(u-/0 ^ G{~u-K), (A.8) 

where u G z[0,ii''], requiring that G{z) has no branch cut on [—6, —a]. The third one comes from the 
assumption that G{z) is regular if |Re(z)| > b. This becomes 

G{u + iK') = G{u-iK'), (A.9) 
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Figure 7: The conformal mapping from z-plane to w-plane. The upper-half z-plane is mapped into the 
region inside the blue rectangle in the u-plane, and the lower-half z-plane is mapped into the region inside 
the red rectangle in the u-plane. 

where u G [—K, K]. 

In |25J, G{u) is extended to the whole u-plane such that the above three equations are satisfied for 
any u G C. It is necessary to check whether such an extension can be done consistently. 

The first step is to extend G{u) to the next region using (jA.Sp . That is, the left-hand side of (jA.Sp 
with u G K + Rq, which is assumed to be already known, define G{u) on through the right-hand 
side. 

The next step is to use 

G{u - 3A") - 2n G{u - K) + G{u + K) = (A.IO) 

which is a consequence of (jA.7l) and (|A.8|) . This equation determines G{u) on Ri provided that G{u) has 
been already defined on and i?;_2j or on and Ri+2- It is easy to check that this extension is 
compatible with (|A.8[) . 

Up to here, G{u) has been defined on IJjg^ 1"^^^ extension to the vertical direction is straightfor- 
ward by using (|A.9p . Obviously, this last step is compatible with (|A.7p and (|A.8[) . 

Now G{u) is defined on the whole u-plane, satisfying (jA.7P(jA.8p()A.9p as functional equations, not 
boundary conditions. The next task is to solve these functional equations. It is convenient to define 

+ ^ ' 2sin(7ri/) ^ ' 

It can be shown that G+{u) satisfies the following equations 

G+iu + 2K) = e-'^*(i-'')G+(u), (A.12) 
G+{u + 2iK') = G+{u). (A.13) 

The original function G{u) can be recovered from G+(w) by 

G(w) = -i [e5'^*''G+(u)-e-^""G+(-u)l . (A.14) 
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One can show that (jA.7|)(IA.8p arc consequences of (|A.12[) while (jA.9p is a consequence of (|A.13[) . There- 
fore, the task is reduced to solving (|A.12[)(|A.13p . 



There exists a solution 

H+{u) 



of ([XT2|)(f03l) where z?q(w) = da (y, (a = 0, 1, 2, 3) are the elliptic theta functions. Note that this 
solution differs from the one given in |25| by the overall sign. One can easily verify that this is a solution 
of (|A.12|(jA.13P provided that e is chosen to be 

£ := i{l - v)K'. (A.16) 

Note that H+{u) has simple poles at u = K and u ^ K + iK' which correspond to z = a and z = b, 
respectively. The conformal map (jA.2[) maps near those points as follows, 



K cx V7^, (A.17) 



u-{K + iK') cx VI^. (A.18) 

Therefore, H^{z) has a branch cut on [a, fo] and it is divergent at the branch points. In addition, iJ+(z) 
has simple zeroes at z = oo (u = iK') and z = e (u = e) where 

e := asn(i(l - i/)A",fc). (A.19) 

In [3S], G+{z) is chosen to be proportional to H+{z). The normalization condition chosen in [35] is 

lim zG{z) = i. (A.20) 

This determines the normahzation constant to be 

= a-^^^H+iz). (A.21) 

This normalization implies the following product formula 

Note that G^{z) has also another branch cut on [— &, —a] but it disappears in G(z). 

The analysis in this paper requires to know the value G.^ (0). Using the identities of elliptic functions, 
one finds 

G+(0) = ie. (A.23) 
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B An identity 



In this Appendix, a6 = 1 is always assumed. 
Consider the fohowing function, 

F{z) := 5(^-i)G+(z-i), ~g{z) := zg{z), (B.l) 

where g{z) is defined in (|4.23p . F{z) has singularities at z = a, 6 just as Gj^{z). One can show that 
g[z^^) has a zero at z = e while G+(z^^) is finite there, implying that F(z) has a zero at z = e. F{z) 
behaves as 

F{z) := —+0{z-^) (B.2) 
ez 

near infinity. 

The product F(z)F{~z) turns out to have a simple form: 

P( ^ _, {z-^-a^){z-'^b^)-^z^^ 

F{z)F(-z) = -z —z IT- G+{z )G+{-z 

(z — e'^y 

_ _, e2(z-2 - a'){z-^ ~ 9) ~ (e^ ^ Q^)(e^ - 9)z-^ 
e2(z-2-e2)(z-2-a2)(z-2-62) 

.,2 „2 

(B.3) 



e2(z2 — a?)[z'^ — 62) 

This implies that z — a,b are the only singularities of F{z), and z — e,oo are the only zeroes of F{z). 
In term of the u variable defined in (IA.2I) . F{z) is written as 

F{u) = g{u + iK')G+{u + iK'). (B.4) 

It is shown in Appendix lAl that G+(u) is defined on the whole u-plane and it satisfies the relations 
(|XT2|)(|XT3l) . Therefore, if g{u) satisfies 

g{u + 2K) = g{u), g{u + 2iK') = g{u) (B.5) 

for any u G C, then, since the ratio F(z)/G-i-{z) turns out to be an elliptic function without poles, it is 
concluded that F{z) is proportional to G'+(z). 
One can show that 

gi-y±iO) - ~giy±iO). (B.6) 
In terms of the u variable, this can be written as 

g{u-K) = g{u + K). {uei[-K',K]) (B.l) 

In addition to this, one has the following condition 

g{u + iK') = g{u - iK'), {u G [-K, K]) (B.8) 

which comes from the regularity outside the branch cuts. These conditions imply that g{u) can be easily 
extended to the whole u-plane while satisfying ()B.5P for any u G C. 
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The proportionality constant is determined from the behavior (|B.2p . This completes the proof of the 
identity 

eg{z-^)G+{z-^) = zG+{z). (B.9) 



C Uniqueness of the solution 



This appendix discusses conditions which uniquely determine the solution of the equation 

/(y) - v{y + iQ)+v{y-iQ)-2nv{-y), (C.l) 
where y G [a, h] C M. The function J{y) can be arbitrary. 



Suppose that there are two solutions vi{z) and V2{z) of (|C.ip . The difference 5v{z) := ^1(2:) — ^2(2) 
satisfies the homogeneous equation 



= 5v{y + iO) + 6v{y - iO) - 2n 6v{-y). 
The general solution of this equation is given in subsection 14.11 The result is 



Sv+{z) = A{z^) + zB{z^)y/{z'^ - a2)(z2 - b^) G+iz) 



(C.2) 



(C.3) 



where 5v^{z) is defined as in (I4.11[) . 
Now, we require that 

• vi{z)^V2{z) are holomorphic on C\[a, &], and 

• vi(z), V2{z) are finite at z = a, b. 

Then, as shown in subsection l4.11 Sv+{z) should be of the form 

dv+{z) ^ \A{z')g{z) + zB{z^)\ ^{z^ - a^){z^ - b^)G+{z), 



(C.4) 



where g{z) is defined as (|4.23p . Both A{z^) and B{z^) are entire functions on C. A further requirement 

• vi(z),V2{z) are finite at infinity 

demands that A{z'^) and B{z'^) must vanish at infinity. Therefore, one finds 

A{z^) = B{z^) = 0. (C.5) 

Since Sv^(z) = implies Sv{z) = 0, it is concluded that the above three conditions determine the solution 
of (jC.ip uniquely. 
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D Determination of coefficients 



The residues of w+(z, ^) at z = ±^ are 



Res^c.+ (z,0 = {c{C)giO+d{m)G+{0, 
Res_5W+(z,0 = {c{Og{-0-d{m)G+{-0, 



where G+{z) is defined as 



The coefficients c(^) and d{^) are the solution of 

gi-OG+{-0 



2z 
d{0 



2(n2 - 1) 



e 2 



62' 



Let the above 2x2 matrix be denoted by A. Its determinant is 

detA = -e(.9(C)+.9(-e))G+(C)G'+(-0 



2C 



and therefore, the inverse is 



2? 



gi-OG+{-0 giOG+iO 



-^G+{~0 -^G+iO 
-gi-OG+i-0 giOG+iO 



The solutions is then 



2(n2 - 


1) 


f 




2(n2 - 


1) 


i 




2(n2 - 


1) 


i 





-e [e-^G+i-O-e^G+iO) 
CGiO- 



2(n2 - 1) e 

^ -r'G(r') 



2(n2 - l)e 



(D.l) 
(D.2) 

(D.3) 
(D.4) 



(D.5) 



(D.6) 



(D.7) 



(Di 
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